Vortex ground state for small arrays of magnetic particles with dipole coupling 



o 

(N 
C 



i 

-f— > 

CO 

c3 



i 

c 

o 
o 



> 
in 

(N 

o 



S. A. Dzian, 1 A. Yu. Galkin, 2 B. A. Ivanov, 3 ' V. E. Kireev, 3 and V. M. Muravyov 3 

1 National Taras Shevchenko University of Kiev, 03127 Kiev, Ukraine 
Institute of Metal Physics, 03142 Kiev, Ukraine 
3 Institute of Magnetism, 03142, Kiev, Ukraine 

We show that a magnetic vortex is the ground state of an array of magnetic particles shaped 
as a hexagonal fragment of a triangular lattice, even for an small number of particles in the array 
N < 100. The vortex core appears and the symmetry of the vortex state changes with the increase of 
the intrinsic magnetic anisotropy of the particle /3; the further increase of f3 leads to the destruction 
of the vortex state. Such vortices can be present in arrays as small in size as dozen of nanometers. 
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I. INTRODUCTION 



Topological defects of vortex type play a paramount 
part in the general physics of ordered media such as su- 
perfluidity, superconductivity and magnetism. In par- 
ticular, vortices and vortex pairs are important in two- 
dimensional (2D) magnetism. Recently, the ground state 
of soft ferromagnetic particles of micron and sub-micron 
size has been shown to be of a vortex type, which has re- 
ceived much attention from the research community^— 
Compared to vortices in supcrfluid systems, magnetic 
vortices in 2D ferromagnets and antiferromagnets have a 
richer behavior, since they may be divided into two dif- 
ferent classes, in-plane and out-of-plane vortices4^ For 
in-plane vortices all spins lie in the vortex plane. Out-of- 
plane vortices have nonzero spin components orthogonal 
to the vortex plane localized within the so-called vortex 
core, a small region near the vortex center. Vortices with 
a core are described by several different types of topo- 
logical charges. 6 " 8 Besides the standard 7Ti-topological 
charge vorticity q, which is similar to circulation in a 
supercurrent systems, one can introduce ^-topological 
charge, the polarity or polarization p — ±1 of the vor- 
tex core, which is the spin direction in the core and is 
connected to the 7r 2 topological charge of the magne- 
tization field. In-plane vortices can be associated with 
the value p = 0. Further, vortices found in the ground 
state of soft magnetic particles should have q = 1 only 
(q = — 1 corresponds to antivortices that can be con- 
nected to "antidots", small holes in a patterned mag- 
netic film 9 ), but they arc additionally classified by the 
discrete number chirality C = ±1, which is the sense 
of rotation of magnetization far from the vortex core. 
For a single vortex in a bulk 2D magnet with easy-plane 
magnetic anisotropy, there is a transition from coreless 
in-plane vortex structure to the vortex with a well de- 
fined core, as the anisotropy strength decreases below a 
certain critical valued The presence of a core plays a 
crucial role in the dynamic properties of magnetic vor- 
tices. In particular, the value of -^-topological charge 
determines special gyroscopic properties of the vortex 
dynamics and the presence of low-frequency dynamics, 
essential for possible applications of magnetic vortices in 



perspective spintronioii — and magnoni o 18 ' 19 devices. 

The presence of vortices in the ground state of soft 
ferromagnetic particles is determined by the balance 
between magnetostatic and exchange energies. For 
disk-shaped particles with negligible magnetocrystalline 
anisotropy and typical thicknesses about 20 4- 50 nm, the 
vortex state is stable if the disk radius exceeds some criti- 
cal value R c ~ 150—200 nm. Vortices in the ground state 
of soft magnetic particles possess a core with the size of 
the order of the exchange length of the material (about 
15 nm in Permalloy). It has been recently shown that 
magnetic vortices of different structure can be the ground 
state for magnetic particles with comparable energies of 
the exchange and dipolar interactions, even for small (of 
the order of 10 3 ) number of magnetic moments in the 
particle. 20 Highly non-uniform ground state is found for 
small particles with a sufficiently high surface anisotropy 

From the perspective of a search for vortex states, the 
promising magnetic systems are those in which the ex- 
change interaction suppressed or non-existent and the 
main source of interaction among structural elements is 
the dipole-dipole interaction of their magnetic moments. 
These are so-called dipolar magnets, i.e. such spin sys- 
tems where a long-range magnetic dipole interaction pre- 
vails. Dipolar magnets have been attracting persistent 
interest during decades as objects of the fundamental 
physics of magnetism possessing some unusual proper- 
ties. One may mention the presence of an ambiguous 
ground state with non-trivial continuous degeneracy even 
for simple cubi o 23 i 24 or 2D square lattices^ - — and the 
existence of special phase transitions induced by an exter- 
nal magnetic fieldi 29 ' 30 Magnon spectra of such systems 
have a non-analytic behavior at small wave vectors^i - — 

The interest in systems with the dominant magnetic 
dipole interaction has significantly increased in last years, 
mainly in the context of artificial magnetic materials 
such as arrays of magnetic nanoparticles. 1 Magnetic sys- 
tems with the dominant dipole interaction possess in- 
teresting physical properties important for applications. 
Among those properties, one can highlight the fact that 
the ground state of an infinite system of magnetic mo- 
ments constituting a lattice and coupled by the dipole- 
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(a) (b) (c) 

FIG. 1: The model of a planar array of magnetic nanoparti- 
cles: (a) the shape of a single elongated magnetic particle, the 
easy axis is denoted as EA, orientation of the magnetic mo- 
ment /2 is denoted by the arrow; (b) a symmetric plane cluster 
comprised of 19 particles; (c) the same in a plane view; mag- 
netic moments of particles in a planar vortex existing at small 
anisotropy are shown with arrows. 



of particles shaped as a fragment of a triangular lattice 
with a high hexagonal symmetry, see Fig. [TJ the ground 
state is a vortex state even for a small number of parti- 
cles in the array. Even for an extremely small array of 
such type, which consists of 7 particles, the vortex state 
energy is almost two times lower than that of a quasi- 
homogeneous state. We show further, that such arrays 
have very interesting behavior: their ground state struc- 
ture is highly sensitive to the change of the anisotropy 
of a single particle. In particular, for a particle made of 
a soft magnetic material this anisotropy could be varied 
by changing the particle shape. 

II. MODEL AND RESULTS 



dipole interaction depends essentially on the lattice struc- 
ture^ and in the presence of the intrinsic (intra-particle) 
anisotropy it depends on the orientation of the easy axis 
of this anisotropy with respect to the lattice axes as well. 

Let us discuss 2D lattices, which will be the object 
of our study. For the lattices of particles with a high 
perpendicular anisotropy, the ground state corresponds 
to various types of two-sublattice antiferromagnetic or- 
der, particularly, a chessboard structure is realized for a 
square lattice^ and a layered one for a triangular lat- 
ticed For finite fragments of such lattices these struc- 
tures vary insignificantly being compared to the infinite 
case. In the case of systems with an in-plane anisotropy, 
the role of the boundaries is more essential. For an infi- 
nite square lattice, the ground state has four-sublattice 
antiferromagnetic order and has a high (continuous) de- 
generacy, while for the triangular lattice a ferromagnetic 
order is realized^ The presence of a boundary, however, 
may change significantly the state of such a system. For 
finite fragments of the square lattice the aforementioned 
continuous degeneracy is removed, but the state remains 
antiferromagnetic ^ 28 ' 30 For a triangular lattice of dipoles, 
the state changes are more radical, and in finite samples 
the ground state can be a vortex state with a closed flux 
of magnetization formed by the magnetic moments ly- 
ing in the plane of the array^ In general, the reason 
of emergence of vortices is the same as for soft magnetic 
dots 3 ^, see also,— ~— where the vortex state emerges due to 
the energy gain of closure of magnetic flux in the sample, 
and the vortices are an alternative to a standard domain 
structure. However, for a rectangular-shaped array the 
vortex state is advantageous only for large enough ar- 
rays. 38 

In view of the importance of vortices for the funda- 
mental physics of magnetism as well as for a variety of 
applications, search for physical systems with vortex-type 
ground state is of great interest. In particular, it is in- 
teresting what is the minimal size of a system carrying a 
vortex, and whether that size can be made significantly 
smaller than the critical size R c indicated above. 

In this paper, we will demonstrate that for an array 



Consider an array of particles placed at the sites of 
a finite hexagonally shaped fragment of a 2D triangular 
lattice. The energy of the this system contains contribu- 
tions from the energy of the magnetic dipole interaction 
and from the energy of the anisotropy: 



where /!,- is the magnetic moment of the particle at site 
I, |/2| = /io, is the magnetic moment of a single par- 
ticle, V = (I — — V\, a is the lattice constant (the 
distance between closest particles in the array plane), 
and (3 is a dimensionless constant determining the mag- 
netic anisotropy strength of a particle. This anisotropy 
is assumed to be uniaxial, of the easy-axis type (so that 
j3 > 0), with the easy axis e z perpendicular to the system 
plane. Here e x ,y,z are unit vectors along coordinate axis. 

One can expect that the absence or presence of a vortex 
core will depend on the effective anisotropy of the system. 
It is important to note that the total magnetic anisotropy 
of the array results from the uniaxial anisotropy of sepa- 
rate particles and from the easy-plane anisotropy of the 
array induced by the demagnetization field of a planar set 
of magnetic moments (similar to the shape anisotropy 
in the thin films). Competition of these two contribu- 
tions determines a complex character of the distribution 
of magnetization in the array. 

If the easy-plane anisotropy induced by the demagne- 
tization field is sufficiently large, that leads to a suppres- 
sion of the vortex core. One can expect that the effective 
anisotropy can be changed by using magnetic particles 
possessing their own intrinsic easy-axis anisotropy. For 
the core to emerge, the effective anisotropy should be re- 
duced, which takes place for particles with the easy axis 
perpendicular to the array plane. For arrays of parti- 
cles made of soft magnetic materials, such situation is 
realized in case of elongated particles oriented perpen- 
dicular to the array plane, see Fig. [T] Such a geometry 
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FIG. 2: Magnetic structure for the cluster of N — 37 particles 
at different values of anisotropy constant f3, here /3i = 0.697, 
02 = 0.729, see the text. The planar components of the mag- 
netic moments for each particle are shown with arrows, and 
vertical moments are represented by circles. For magnetic 
moments with the positive or zero z-projection, open sym- 
bols (arrows with open heads or open circles) are used, while 
particles with negative z-projection of magnetic moments are 
depicted by the solid symbols. 



of the problem is promising for ultradense information 
storage^— and is naturally realized, for example, if the 
array is created by self-organization of small elongated 
particles floating in liquid. Competition of these two 
magnetic interactions provides the possibility to change 
the effective anisotropy of the system, which, as we will 
show below, allows one to impact the structure of the 
magnetic macrovortex. 

To analyze this problem, we have employed two meth- 
ods: numerical minimization of the energy ([1]) using 
the standard Gauss-Seidel algorithm, as in^S and Monte 
Carlo analysis using the simulated annealing technique. 43 
The energy minimization has been performed as follows: 
we start with j3 = 0, choose a simple in-plane vortex as 
an initial condition, and numerically minimize the en- 
ergy then the value of (3 is increased step by step. This 
method is working pretty fast and gives a good descrip- 
tion of the structure for continuous transitions, see be- 
low, while the Monte Carlo method is important for the 
analysis of points where the magnetic structure of the 
system is changing discontinuously, with a coexistence of 
(metastable) states close to the transition point. 

Numerical calculations have been carried out for com- 
paratively small clusters shaped in the form of a regular 
hexagonal fragment of the lattice consisting of 19, 37, 
61, 91, and 127 nanoparticles. For all studied systems, 
we have found an in-plane vortex in the ground state at 
small enough /3. The structure of the vortex changes con- 
siderably as (3 increases and passes through two critical 
values /3i and #2, see. Fig. [2] and the detailed discussion 
below. Further, we have found a prominent transition at 
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FIG. 3: The dependence of the energy W on the anisotropy 
constant j3 for clusters of different sizes (labels in the legend 
denote the number of spins in the cluster), found by numerical 
minimization of ([l]). Horizontal lines denote the energy of 
three-domain antiferromagnetic state, see Fig. |2(e)"| while the 
sloped curves correspond to the energy of a vortex state. 



some value (3 = /3 t > from the vortex state to the 
state with a fragment of the antiferromagnetic structure, 
and close to f3 — f3 t we have observed a noticeable region 
of coexistence of vortex and antiferromagnetic states, see 
Fig. [3] The behavior of the energy as a function of the 
anisotropy constant j3 near this transition is similar to 
that for a thermodynamic potential as a function of tem- 
perature near first-order phase transition. On the other 
hand, the dependence of the energy on f3 did not exhibit 
visible peculiarities at j3 = /3i,2, where a change of the 
vortex structure has been detected. 

To analyze the structure and symmetry of the vortex 
core, after completing the energy minimization for each 
value of j3, we have been calculating the value of the out- 
of plane component M z of the total magnetic moment, 
as well as the length of the planar component M p \ = 



+ My. It turns out that just these parameters are 

most sensitive to the vortex structure and allow observing 
peculiarities of the vortex core behavior, see Fig. |U 

The behavior of the total magnetic moment of a par- 
ticle array with a vortex is rather complicated. At small 
anisotropy, the picture remains the same as for isotropic 
particles^ and a vortex with purely planar distribution 
of magnetic moments is realized. In this case, accord- 
ingly, the total z-projection of the magnetic moment van- 
ishes, but the planar component of the total magnetic 
moment M p \ is non-zero. With increasing f3, a nonzero 
value of M z emerges at some critical anisotropy j3 = fi\. 
The fact that M z ^ means that an out-of-plane core ap- 
pears. However, with the appearance of nonzero M z , the 
planar component of the total moment does not vanish 
immediately, i.e. within some finite interval /3± < ft < P2 
both in-plane and out-of-plane components of the mo- 
ment are nonzero, see Fig. 2] The planar component van- 
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FIG. 4: The in-plane and out-of-plane magnetic moment com- 
ponents M z (full symbols) and M p i (open symbols) vs the 
anisotropy constant ft, in units of po, for a cluster with 61 
particles. The details of behavior in the region of vortex core 
reconstruction, fti < ft 2 < ft, here fti ~ 0.744 /? 2 ~ 0.786, are 
presented in Fig. [5] on a different scale. Up triangles label the 
value M z — no for the antiferromagnetic state. 

ishes for ft > ft 2 , whereupon the vortex state structure 
becomes more symmetric than that observed at ft < ft 2 , 
see Fig. [5] Such a symmetric vortex structure is ob- 
served within a wide range of ft, ft 2 < ft < ft t and the 
magnetic moment M z changes considerably while M p i re- 
mains zero. With the further increase of the anisotropy, 
the vortex structure becomes an antiferromagnetic struc- 
ture similar to that which is found for strong perpendic- 
ular anisotropy. The presence of boundaries leads to the 
emergence of three different domains of such structure, 
as dictated by the system symmetry. 

Note that the behavior of M z and M v \ is similar to that 
of order parameters near second-order phase transitions. 
This observation can be used to perform a symmetry 
analysis of the transitions between the different vortex 
states. It is important to obtain analytical results, due 
to the limited accuracy of numerical data, but also be- 
cause the numerical analysis is hindered near transition 
points ft — fti.2 because of the "critical slowing down" 
of relaxation similar to that found near a second-order 
phase transition, which manifests itself in a substantial 
increase in the numerical calculation time. Therefore, we 
should study the possibilities of existence of such transi- 
tions in our system from the viewpoint of symmetry. 

III. SYMMETRY ANALYSIS 

To describe the complex character of changing the vor- 
tex core structure let us use symmetry arguments in line 
with the phase transition theory of Landau. For both 
observed critical values of anisotropy, the symmetry of 
state changes essentially. At ft > fti, when out-of-plane 



FIG. 5: The in-plane and out-of-plane magnetic moment com- 
ponents M z (up triangles) and M p \ (down triangles) vs the 
anisotropy constant ft, in units of /io, for a cluster with 61 par- 
ticles, in the region of the vortex core reconstruction. Note 
that the scales for M p i and M z differ by the factor of 100. 



core emerges for the first time, the sign of M z can be ar- 
bitrary, i.e., there is a spontaneous breaking of Z2 sym- 
metry with respect to M z at ft = ft\. In contrast to that, 
the symmetry of the planar distribution of magnetic mo- 
ments does not change at this transition point: within 
the range ft\ < ft < ft 2 it remains the same as for ft < fti . 
At the other transition point ft = ft2 the situation is 
different: the dependence of the out-of-plane component 
M z (ft) does not have any visible peculiarities, while the 
planar component M p \ vanishes at ft = fti and remains 
zero for ft 2 < ft < ftt, i.e., up to the point of destruction 
of the vortex state at ft — ftt, see the detailed graph in 
Fig. El 

For the ground state in the interval ft 2 < ft < ftt, we 
observe a higher symmetry of the moment distribution 
than the vortex states with M p i ^ at ft < ft 2 (specif- 
ically, within the numerical accuracy of our simulations, 
we observe the C'e symmetry; it would be worth finding 
out whether this symmetry is exact). Rather low sym- 
metry of the vortex ground state at ft < ft 2 is caused by 
the presence of nonzero M p \, which can be traced down 
to the presence of a non-zero planar component of the 
central magnetic moment fto- Obviously, the emergence 
of non-zero M z at ft > fti is caused by fa coming out of 
plane. Writing down /Io as fto = /io (sin6*oe 2 + cos6>oe p i), 
where e p \ lies in plane of the system, we obtain that the 
transition at fti is connected to the appearance of non- 
zero 6> , with 8 = at ft < fti and 6 ^ at ft > fti. As 
we pointed out, the symmetry of the state with M z 7^ 
is lower than for a planar vortex, therefore the value of 8q 
serves as the order parameter for the transition at ft — fti . 
In this case, one can expect that at ft > fti the behavior 
of the "order parameter" close to the transition is given 
by #0 oc v 'ft — fti and is characterized by singular behav- 
ior, ddo/dft — > 00 at ft — > fti + 0. On the other hand, the 
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FIG. 6: The value of the angle 9o (normalized by it/2) vs. 
the anisotropy constant j3 for a cluster with 61 particles in 
the region of the vortex core reconstruction. Symbols present 
the results of numerical calculations, while the dashed and 
solid lines show analytical results from the phenomenological 
theory. The dashed line corresponds to Eq. ([3]), and the solid 
line shows the result that follows from including corrections 
due to the sixth-order term with /3, see the text. 



hexagonal symmetry in the spin distribution may appear 
only when the central moment is directed strictly per- 
pendicular to the system plane, i.e. at 9q = n/2. If at 
P = 02 the symmetry increases up to the hexagonal one, 
the quantity $o = tt/2 — da should serve as the order pa- 
rameter for this transition. We arrive at the conclusion 
that the behavior of the out-of-plane component of the 
central spin /io, i.e., the dependence #o(/3), dictates the 
change of symmetry of the vortex state. The information 
about the full 9q(P) dependence can be obtained only nu- 
merically, but the presence of square-root singularities at 
P — > fii + and j3 — > p 2 — is rather easily verified, see 
Fig.! 

The detailed analysis of the 9q(P) dependence allowes 
one to present a closed phenomenological expression for 
the "thermodynamic potential" $ that defines the be- 
havior of 0q{(3) in a wide range of p. Indeed, in line with 
the Landau theory, this potential can be constructed in 
the form of the expansion in powers of the order param- 
eters, which are 9 &t ft ~ fli or z9 = ir/2 — 8 at (i ~ /3 2 - 
Equivalently, sin 9q and cos 9q can be used instead of an- 
gles 9 and d . As odd degrees of [i 0z = sm8 are for- 
bidden by the condition of the time reversal invariance, 
and the simplest form of this energy is the following: 
$ = A sin 2 6q + B sin 4 9q- It is easy to see that, up to an 
inessential overall factor, the correct behavior is provided 
by the expression 



$ = 2 (A - P) sin2 d ° + 4 (A " A) sin4 °> ( 2 ) 



FIG. 7: The dependence of /?2 (up triangles) and /3i (down tri- 
angles) on the number of particles in the cluster, lines present 
the dependence fitted by the logarithmic function, see the 
text. 



which leads to the simple result: 



an0 Q = J£—BL, Pi<p< fa, (3) 
y P2 — Pi 

while 6» = at f3 < f) u and 9 = vr/2 at f3 > /3 2 . 

Such simple dependence describes the numerical data 
fairly well, see Fig. [BJ Deviation from the simple law 
given by (J3)) can be accounted for by adding the term 
{(3/6) sin 6 6*o to the expansion ([2]). As seen from Fig. El 
it provides a perfect description of the numerical data at 
sufficiently small f3, typical values are j3 < 0.1 (02 — Pi)- 

The critical values of the anisotropy constant /3i and 0% 
grow with the increase of the cluster size N, see Fig. [71 
though for the studied values of N this dependence is 
rather slow. The numerical data for N > 37 can be well 
fitted by a logarithmic dependence of the form = 
Ai, 2 + B h2 \nN, where A x = 0.38397, B x = 0.08703; 
A 2 = 0.36236, B 2 = 0.10222. 

The value of #o not only dictates the vortex core sym- 
metry but it also quantitatively defines the important 
vortex characteristic, the total out-of-plane moment of 
the particle with the vortex. Singularities in the #o b e_ 
havior at ft = pi are reflected in the M z (P) depen- 
dence, M z oc sin 9q near this point. It is worth noting 
that M z plays a special role in the dynamic properties, 
namely, M z serves as a proper collective variable describ- 
ing the radial mode (with the azimuthal number m = 0) 
of the magnetization oscillations in the particle with a 
vortexj^i the theory is in agreement with recent experi- 
menti 4 ^ Therefore, the presence of singularities in M Z (P) 
should manifest itself in the behavior of an equivalent of 
this mode for the considered system. 

In addition, the presence or absence of the vortex core 
is also important for the properties of azimuthal modes 
with m = zfcl i 46 ! 47 For a purely planar vortex, the modes 
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with m = ±1 form degenerate doublets, while the emer- 
gence of the core leads to splitting of these doublets. 
Thus, one can expect a crucial impact of the vortex struc- 
ture modification on the properties of eigenmodes of the 
vortex-state particle, although a detailed discussion of 
the dynamical properties of the particle with the vortex 
is beyond the scope of this work. 

IV. SUMMARY AND DISCUSSION 

To conclude, we have shown that high-symmetry 
hexagonal fragments of a 2D closely-packed triangular 
lattice of magnetic particles contain a vortex in the 
ground state, even for a small fragment size. The vor- 
tex structure is very sensitive to the intrinsic anisotropy 
/3 of the particle. At small anisotropy, there is a purely 
planar vortex. With the increase of j3, the symmetry of 
the vortex ground state lowers initially at some critical 
value /3 = 01, and then increases to a high sixfold axial 
symmetry above another critical value /3 = f$i > f3\. It is 
worth to note that those two transformations bear a simi- 
larity to second-order phase transitions. Both transitions 
take place at sufficiently weak anisotropy, the dimension- 
less parameters /3i.2 do not exceed one. This value is 
essentially smaller than the easy-plane anisotropy of a 
planar array induced by the demagnetization field with 
the characteristic value /3 a rray ~ 10, see Ref. Ac- 
tually, this anisotropy is smaller than it is necessary to 
create the perpendicular magnetization of a cylindrical 
magnetic dot. 

An important challenge in the physics of magnetic vor- 
tices is to find ultra-small (smaller than 100 nm) systems 
with vortices in the ground state, this problem is of great 
interest for both fundamental physics and applications. 
In addition to lithographic magnetic materials, where the 
particle size is of the order of tens of nanometers ^2r— the 
proposed theory is applicable to other 2D systems with 
anisotropic particles having magnetic or electrical dipole 
moment. - 7 One could expect that if an array can be com- 
posed from small enough particles, having finite magnetic 
or electric dipole moment, the vortex state will be present 
for arrays 10-20 times larger than the particle size. Such 
systems can be realized for composite magnetic materi- 
als, for example, for granular magnets with the content of 
the magnetic component less than the percolation thresh- 



old, where the exchange interaction between nanometer- 
sized grains is anomalously small. Another example is 
the inhomogeneous state arising in the vicinity of the 
metal-insulator transition in doped manganites, which 
involves small particles of the ferromagnetic (metallic) 
phase distributed over a nonmagnetic host; their phys- 
ical properties are determined to a large extent by the 
dipolar interactions between these particles 4£ The ex- 
perimental implementation of the artificial crystals, in 
which particles with magnetic moments of the order of 
10 3 Bohr magnetons form an ordered lattices, has been 
reported recently^ As one more example, it is instruc- 
tive to mention a new class of materials, namely, molec- 
ular crystals formed by high-spin molecules. The total 
magnetic moment of such a molecule can be as high as 
dozens of Bohr magnetons, but the exchange interac- 
tion between magnetic moments of different molecules 
is almost negligible^ Note also so-called dense phases 
formed by nanometer-sized magnetic particles moving 
freely in a liquid (that is the standard situation for a 
ferrofluids). 51 For all these systems with a particle size 
of the order of nanometers the vortices described here 
can be present for objects as small as dozen of nanome- 
ters; those are, to the best of our knowledge, the smallest 
vortex-bearing systems discussed in the literature. 

It is worth noting that the presence of a vortex ground 
state for such small systems and the transitions with the 
vortex core reconstruction is a consequence of the high 
(hexagonal) symmetry of the array. For square or rect- 
angular arrays the vortex state appears for large enough 
arrays only, 3 - As we found, for an array shaped as a 
regular triangle, the vortex state could be present for 
small arrays, but with the increase of the anisotropy 
the vortex remain coreless all the way till the transi- 
tion to antiferromagnetic state. The two-dimensional 
nature of an array is also quite important. Thus, two- 
dimensional closely-packed arrays of magnetic particles 
represent vortex-bearing systems with potentially small 
sizes and offer a unique possibility for manipulating the 
symmetry and structure of the vortex core. 

We are thankful to V. G. Baryakhtar, A. K. Kolczhuk 
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